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Abstract 

Non-equilibrium modifications of the nuclear equation of state are assessed 
"q ■ in the framework of Extended Irreversible Thermodynamics. The expected 

size of the non-equilibrium corrections to quantities like the nuclear-matter 
incompressibility K and the liquid-vapor critical temperature T c are com- 

X 

pared to the uncertainties associated at present with the usual equilibrium 
values. The quantitative estimates are obtained in situations characterized 
by a stationary planar shear flow of realistic magnitude, or by a stationary 
temperature gradient. 
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The study of nuclear collisions became a major subfield of nuclear physics in the past 
two decades. Due to the complicated nature of the physics and due to the many degrees of 
freedom involved, tools from statistical physics and thermodynamics proved to be efficient 
in modeling and analyzing these collisions. For example, the concept of nuclear temperature 
is commonly used. However, the applicability of equilibrium concepts like temperature to 
the nuclear collision scenario is frequently questioned. 

It is hard to believe that equilibrium is reached during the short lifetime of the collision 
system in a high-energy nuclear collision, and thus the above question is well-justified. One 
answer, which tries to avoid the problem of equilibration in each collision, points out that 
these experiments deal with very large samples of identically prepared collision systems. 
What we have then, is an ensemble of events, while the usual introduction to statistical 
physics proceeds via postulating the identity of time averages with ensemble averages (er- 
godicity). In nuclear collisions, we deal with these ensembles directly. Each member of the 
ensemble need not be in equilibrium, and the quantities used to characterize the ensemble 
(e.g. temperature) are ensemble averages. Due to the relatively small number of degrees 
of freedom (compared to e.g. Avogadro's number), the distributions are in fact much wider 
than in conventional statistical physics, leading to the important role of fluctuations in these 
collisions. 

Many people find the above argument unsatisfactory. In this Letter we consider an alter- 
native point of view. We want to be able to use statistical-physics concepts like temperature 
without relying on the assumption of equilibrium. One may distinguish between global equi- 
librium of the system (or a well-defined subsystem) and local equilibrium (such as in the 
viscous hydrodynamic model). Here we are particularly interested in the situation when not 
even local equilibrium is established. The question then becomes: is it possible to describe 
and usefully quantify the deviation from equilibrium? Similar problems are addressed in 
non-equilibrium thermodynamics In particular, attempts to study non-equilibrium flu- 
ids have been made in the framework of "Extended Irreversible Thermodynamics" (EIT) 

EH- 
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EIT starts with the assumption of a local thermodynamic potential (see below), which 
depends (in addition to its usual variables) on the dissipative fluxes in the system such as 
the heat flux or the viscous part of the pressure (or momentum-flux) tensor. The dissipative 
fluxes are taken as additional independent variables. While the evolution of the standard 
variables is governed by the conservation equations of mass, energy and momentum, EIT 
derives further equations for the dissipative fluxes. These equations can be used to calculate 
non-equilibrium corrections to such quantities as the isothermal nuclear-matter incompress- 
ibility, K = 9/(pKx), where is the compressibility usually defined as kt = p 1 (dp/dp)T- 
Although equilibrium in the common sense is not assumed within EIT, systems are easiest 
to discuss when stationary. The examples we wish to consider are those of a planar flow 
under shear, and of a constant temperature gradient. The nuclear collision scenario is of 
course far from stationary. We believe, though, that an estimate under these assumptions 
will indicate the order of magnitude of the necessary non-equilibrium corrections. 

The value of K for nuclear matter has been subject to intense debate. Values have been 
extracted from the giant monopole resonance supernova calculations ||, and from 

nuclear collisions assuming momentum-independent ]7||| and momentum- dependent [PHTT|] 
mean fields. In general, the incompressibility (isothermal or adiabatic) is expected to be 
a function of density and temperature. Numerical values quoted in the nuclear context 
normally refer to standard density and zero temperature. Recently, information from several 
independent experiments on sideward flow was assessed within the framework of a transport 



model to arrive at 165 MeV < K < 220 MeV fll2| . We will approximately represent this 
as K = 190 ± 30 MeV (see also |13[)- We assume that ±30 MeV is the characteristic 
uncertainty associated with the equilibrium value of K in a reasonable range of densities 
and temperatures. While the value of K may be further discussed on the basis of different 
calculations, to try to achieve better precision in similar studies appears to require a very 
large effort. 

In what follows we focus on the uncertainty in the value of the isothermal incompressibil- 
ity. Note that the extracted values of K refer to nuclear matter in equilibrium, even though 



transport models do not rely on equilibrium. It is interesting to compare the estimated 
EIT non-equilibrium correction to the quoted uncertainty: significant corrections would in- 
dicate that the equilibrium value of K is not really relevant in nuclear collisions. If the 
non-equilibrium corrections are comparable to the uncertainty, then efforts to determine K 
with better precision are not likely to be meaningful. 

In order to estimate the non-equilibrium contributions to the incompressibility or other 
quantities, specific situations must be considered. Within EIT PH, the complete differential 
of the Helmholtz free energy per particle is represented locally as 



df 



-s dT — pdp 1 + a q ■ dq + j3 Ii v : dU v 



(1) 



where U v is a traceless viscous pressure tensor proportional to viscosity 77 and to velocity 
derivatives, double dot indicates double contraction, q = kVT is the heat flux, and k stands 
for heat conduction. As the nuclear incompressibility can be determined from sideward flow 
in energetic (few hundred MeV/nucleon) collisions, where shear develops at the center of 
a colliding system H], we choose the situation of a planar shear, characterized by 
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with V = i] (du/dx), where u is the velocity along the y axis (x and y axes perpendicular 
to the beam direction z), for assessing the corresponding correction. The coefficient in the 
shear term in (P is estimated as 



2r] P 



(3) 



where r v is the relaxation time for n„. For a rarefied gas, with t v ~ (pcv) 1 and 77 ~ mv/3a, 
the coefficient reduces to 



P 



1 

2p 2 T ' 



(4) 
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see also 0. From the equality of mixed derivatives of the Helmholtz energy, assuming 
an inverse quadratic dependence of f3 on density, the pressure under shear alone becomes 

p = p -(3p% :% = Po -2(3pV 2 , (5) 

where po = po(p, T) is the pressure for given p and T in the absence of shear. 
Using (HI) we find for the incompressibility (see also Hl5| ) 



A '= 9 (l)^ A '° +i8/3p2 ' (e) 

where Kq = 9(dpo/dp)T is the equilibrium value. Thus the incompressibility increases under 
shear. The larger is the density, cf. (f|), for a given 7 = du/dx, termed shear rate, the lower 
is the nonequilibrium correction, consistent with the shortening of the mean free path. 

We next assess the actual magnitude of the correction to the incompressibility under 
typical circumstances in collisions. For this purpose we use results of a recent transport 



simulation. Figure 11 of Ref. |n| shows density, velocity-component, and entropy profiles 
across the center of a Au + Au system at 400 MeV/nucleon, at an intermediate impact 
parameter, during the compression stage. The shear rate at the center of the reaction in 
the figure is 7 » 0.07c and the density is p m 0.30 fm -3 . Given the density and entropy 
there, we estimate a temperature of T « 44 MeV. The nuclear viscosity for these parameter 
values is about 77 ~ 55 MeV/fm 2 c f!6 |. Combining the above results and Eqs. (|4]) and (||), 
we obtain K — Kq ~ 20 MeV. Similar corrections are obtained at lower beam energies, 
as long as effects of Fermi degeneracy do not yet enter the picture. We conclude that 
the nonequilibrium correction is comparable to the present uncertainty associated with the 
determinations of the incompressibility. 

Another quantity of interest, the critical temperature of the nuclear liquid- vapor phase 
transition p7}-|T9|j can be used as a further illustration of the corrections to the equilib- 



rium values obtained in EIT. This transition is explored in low-energy reactions, in central 
and/or in spectator regions, where shear may be low, but sizeable temperature differences 
are possible between the matter exposed to and shielded from vacuum. The outward heat 



flow can modify local thermodynamic quantities and equations of state in a manner similar 
to the case of shear. Both shear and heat flow are known to reduce critical temperatures 
in macroscopic media |[20|| . For this illustration, we adopt a simple equation of state of the 
Van der Waals type for nuclear matter in the transition region, and consider a constant 
temperature gradient. 

The coefficient in the heat term in (]l|) may be estimated as [20 



where r q is the relaxation time for q, and, for a rarefied gas with k ~ c p v/3a, 

2 m 

see also 0. Under an assumption of inverse quadratic dependence of a on density, the 
pressure, for a given heat flux, becomes 

p = p -apq 2 i (9) 

where p = p (p,T) is the equilibrium pressure. Insertion of this result into the Van der 
Waals equation of state then yields 

(p+y 2 + ^f) {y-b)=T, (io) 

where v = 1/p is the specific volume. The critical-point conditions 



give then a lowering of the critical temperature p0 |. To leading order in q 2 , the tempera- 
ture is 

T c {q)=T° c -\aq 2 . (12) 

The phase-transition region is crossed in reactions when matter expands to the vacuum 
and nuclear fragments are formed. If the expansion is isentropic and statistical effects 
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are ignored, then the product pT -3 / 2 may be considered constant. Further, if the density 
changes from p ~ p i n the interior of the expanding matter to p ~ 0.2 p at an outside freeze- 
out surface, over a characteristic distance ~ (5 — 10) fm, then the resulting temperature 
gradients are of the order of 1 MeV/fm in the critical region T c ~ 15 MeV. Given that 
actual gradients would be somewhat moderated by the heat conduction (as is supported 
by numerical simulations JT^), we expect heat fluxes ]16[ of magnitude g~0.1 MeVc/fm 3 . 
With (|8|) and p ~ 0.5 po in the critical region, we finally estimate the correction to the 
critical temperature to be of the order of ~ 2 MeV. 



On examining such corrections as in ( |T2"D or (|5|) for a rarefied gas, it is seen that these 
corrections are associated with the variation of temperature or of collective velocity over 
the distance of one mean free path. It is plausible that the critical temperature of un- 
charged infinite matter could be less relevant for fragment formation, than e.g. the limiting 



temperature pi| , |22| , the maximum temperature a residual nucleus could withstand given 
Coulomb interactions and the lowering of surface tension. In any case, EIT gives a typical 
magnitude of corrections on account of dissipative fluxes, and these corrections are compa- 
rable to typical uncertainties associated with the choice of interactions [p3fl . 

To summarize, in this Letter we quantified the deviations from equilibrium expected 
in intermediate energy nuclear collisions. To achieve this goal, we used the tools of Ex- 
tended Irreversible Thermodynamics. As specific examples, the non-equilibrium corrections 
to the incompressibility and to the critical temperature of the liquid-vapor phase transition 
were examined. Our numerical estimates are based on a stationary planar shear flow and 
on a constant temperature gradient, respectively. The comparison of the non-equilibrium 
modifications to the uncertainties associated with equilibrium values reveals that these are 
comparable to each other. This casts doubt on the reducibility of the uncertainties on the 
basis of nuclear collision data. The nonequilibrium modifications imply that nuclear matter 
acts effectively stiffer and displays a lower critical temperature than in equilibrium. 
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